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The  Aerobaliistio  Re-entry  Technology  (ART)  program  of 
the  Haval  Surface  Weapons  Center  has  an  objective  to  provide 
the  aeroballistic  tecnnology  required  for  the  design  and 
development  of  advanced  re-entry  body  concepts.  Essential  to 
such  an  effort  is  the  ability  to  predict,  evaluate  and  modify 
aerodynamic  characteristics  of  candidate  con? iguraticns . Ground 
and  flight  test  procedures,  as  sell  as  computer  modeling  and 
analysis,  are  used  in  the  general  approach.  Full-scale  flight 


testing  is  compile at 
testing  as  currently 


y ana  insecure. 


T aflnl  n: 


most  practical  way  to  oi 


•oaynaiaac 


data.  The  ballistic  range  ahd  the  hypervelocity  wind  tunnel 
are  the  most  frequently  used  of  the  ground  facilities  in  BV 
data  gathering.  Ballistic  ranges  have  the  obvious  advantage  of 
permitting  measurements  to  be  made  in  the  absence  of  supports 
and  their  associated  interference  cf  the  flow  about  the  test 
body.  However,  range  testing  is  costly  on  a per  shot  basis, 
flight  parameters  are  difficult  to  control  and  data  acquisition 
is  very  complicated.  Nevertheless , it  has  been  suggested 
that  in  wind-tunnel  testing;  the  flow  interference  due  to  supports 
results  in  uncertainties  in  /the  final  results  . 1 In  order  to 


determine  if  the  Hyperveloeity  .Facilities  at  NAVSUKFWFHCEH  are 
usable  tools  for  studying  re-entry  body  dynamic  stability, 
a test  program  was  undertaken . 

The  Supersonic*  Tunnel  Association  standard  cone  with, 
varying  nose  bluntness*  was  selected  for  testing.  This 
configuration,  a 10-degree  half-angle  done,  was  chosen  because 
it  Is  representative  of  a classical  re-entry  .shape  and  because- 
extensive  data  exist  fhf  this  conf igurat  1 on -at  lower  Ms  ch 
numbers.  However,  no  dynamic  Sat*  are  available  for  this  _/ 


stanoara  configuration  at 


jerquing  wacn  number  tli  = ic; 


w > «*»•?  WU.»*  W*  O W W»  v Vi*  i M4*<*w*  ■ 

Fitch  damping  measurements  were  made  as  a function  of 
angle-Gf-attaek,  nose  bluntness,  Reynolds  number  and  Mach 
number  using  the  forced  and  free,  oscillation  techniques.  Data 
thus  obtained  are  then  compared  with  data  from  other  facilities 


"Welsh,  G.  J.,  Winchenbach , G.  L.  and  Hadagan,  A.  , 
"Free  Flight  Investigation  of'  the  Aerodynamic  Characteristics 
of  a Cone  at  High  Mach  Hashers,"  AIA.A  Journal,  Vol.  8, 

Ho.  2,  February  1970,  op.  29^-300 
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ivj  ~.ased  upon  differences  between  In-plane  and  out- 


While  good  agreement  Is  noted  with  other  wind  tunnel  data, 
there  appears  to  be  a signif leant  discrepancy  between  these 
wind-tunnel  measurements  and  the  trend;  Indicated  by  the 

range  data.  This  paper  will  offer  an  explanation  for  this 

-? 

• iV  J u*/Wii  u 1 ' * 

of- plane 

rsnrS  ts 

he  presentation  format  will 

a derivation  of  the  equation  of  motion  of  a one-degree- 

of- freedom  system  to 


n * i myyg 


inc  1 


*.  Out-of-plane  damping  is  present  in  the 
absent  in  the  single-degree-of-freedom 
iurements . 


■ble  analysis  of 


negan  of  the  Naval 


>ons  Center  are  heavilv  relied  us on . 


Nest , the  data  reduction  technique  Is  discussed  including 
a description  of  the  Differential  Correction  Procedure  and  a 

Fortrav  XV  computer  program  utilised.  Here  the 


H HIT 

work  of 


Ikenherr/'*  and  Nielsen5  have  contributed. 


The  one-degree-of- freedom  free  and  forced  oscillation 
systems  are  described  and  the  wind  tunnel  tests  conducted  are 
outlined  and  results  presented. 

Finally,  ah  explanation  of  the  range  and  -wind  tunnel 
damping  discrepancy  is  presented.  This  explanation  is  based 
on  work  recently  completed  at  the  Aerospace  Research  Laboratory. 


“Orlik-HucKemann,  0.,  "Wind  Tunnel  Measurements  of  Dynamic 
Derivatives j " lecture  Motes,  University  of  Los  Angeles,  August 
(Available  from  National  Research  Council  of  Canada,  Ottowa) 
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DEVELOPMENT 


EQUATIONS  OF  MOTION 


The  equation  of  notion  of  a viscously  damped,  one-decree- 
of-freedcm  oscillatory  system  -an  be  expressed  as:  rr**®|**^' 


Ik  U) 


for  the  forced  oscillation  case.  If  no  external  forcing 
functions  are  present;,  the  system  takes  the  fora  of  a free 
oscillation  represented  by  the  Homogenous  equation: 

Mx+Ck+  K*  m Q 

Tn  the  above  expression: 

M is  the  mass  Kx  is  the  inertial  reaction 

C is  the  damping  coefficient  Cx  is  the  damping  reaction 

Sis  the  spring  coefficient-  Kx  is  the  spring  reaction 

? is  the  forcing  function  amplitude  Bcoswt  is  the  forcing  funetio 

x:  is  the  generalized  coordinate  « is  the  frequency  of  the 

forcing  function” 


Consider  a wind  tunnel  moqei  angle-of-attack  (see  Figure  30s, 
If  the  model  was  not  supported;  yet  constrained  to  one-Cegrees 
ofrfreedome  motion  a pitchifis  moment,  due  to  the  pressure 
induced  normal  force  and  a diaplng  moment,  always  opposing’ 
oscillation  would  be  generated.  If  the  model  were  supported;, 
such  oscillations  would  occur  about  a trim  angle-of-attack- 
This  trim  angle-of-attack  is  representative  of  a balance 
between  the  pitching  moment  and  the  moment  generated  by  the 
support  in  response  to  its  angular  displacement  (see 
Figure  2 ). 

The  instantaneous  amplitude  of  the  model  about  the  trim 
angle,,  a-,.  Is  represented  by: 

«i.s  et~  «<_  (3) 


For  this  trim  angle  case,  equation  (1)  can  be  used  to  model 
the  interaction  of  the  moments  actinE. 


MOMENTS  ACTING  ON  WIND  TUNNEL  MODEL  AT  ANGLE-OF-ATTACK 
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Where  Iv, 


, M and  M0  are  the  moment  of  inertia,  dynamic 
"peed,  amplitude  of  time-dependent  applied  moment 


y>  '** 

pressure^  air 

and  trim  moment  required  to  maintain  the  trim  angle.  Consider 
the  condition  of  o<<  1.  Under  this  assumption,  we  can  expand 


the  damping  derivative,  Cmq,  and  the  pitching  moment  term,  C^, 


in  a series  about  the  trim  angle,  ctm^  4 _ 

liWLrl 


= M SiH&4%+  Mb  6 *¥•** 

Considering  equation  (^)  for  the  case  of  steady  state  oscillation: 


(5) 


-9*accM(<>l * m6 


<>t 


G 


(6) 


*V 

* 


y t 


(7) 


which  is  representative  of  the  pitch  moment,  support  displacement 
generated  moment  previously  mentioned.  As  the  damping  moment 
is  an  even  function  of  a,  its  odd  derivatives  must  vanish.  The 
pitching  moment  is  an  odd  function  of  q,  thus  its  even  derivatives 
must  vanish,  including  the  above-mentioned  Considerations, 
equation  (5)  necomes:  „ _ _ 

. a.  & mat 

^ *<r  g-  - g^-ljjsg^  _ . 

If  the  angle  of  oscillation,  a,  is  sufficiently  small,  equation  (7) 
becomes : 

x d ^ 

Sin  tt>4t 


#• 


•it. 


•<  s 


In  order  to  be  more  concise,  the  equation  (.8)  may  be  written 
in  the  form: 


d*5  rsaiV  «VWc  la  - 

^ "T*  * - t * - r StHiOfX. 


(9) 
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The  damping  term.  Mil  > contains  both  aerodynamic  and  structural- 
mechanical  damping  and  the  stiffness  term, Mi  > contains  both 
aerodynamic  and  mechanical  stiffness. 

The  homogeneous  form  of  equation  (9 ) may  be  solved  using 
the  Laplace  transform € is: 

m - * 


with  initial  conditions: 


alt..-* 

Equation  (9)  then  becomes: 


of  rearranging  yields: 

<*.  is 


So  cs- 


< 

_ s i 

raking  the  inverse  transform  of  equation  (12)  yields: 


oscillatory  motion  ceases);  when M^gfothe motion  is  undamped  and 
the  motion  oscillates  at  a frequency  of  J* w , These  two 


extremes  of  critically  damped  and  undamped  motion  may  be 
represented  by  the  following  parameters: 


- m/x 


6 Mixon,  F.  E. 
Englewood  Cli f fs , 


, Handbook  of  Laplace  Transforms, 
New  Jersey,  19^0 


Prentice-Hal 


(10) 

(ID 

(12) 

(13) 

(D) 

(15) 

l,  Inc., 
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The  damping  ratio  represents  the  amount  of  damping  present 
aa  a decimal  fraction  of  the  amount  required  for  critical 
damping.  . t - — , 

= /-MW  I 


(16) 


is  the  undamped  natural  frequency. 


Thus  equations  (13) 


and  (14)  may  be  rewritten: 

siHtTTVW 


(17) 


(18) 


Equations  (10)  and  (11)  can  also  be  used  to  rewrite  equation  (9)  as 


= P.?m 


2 . 


(19) 


Consider  equation  (17).  The  amplitude  of  the  t th  oscillation 
can  be  represented  as: 

• JhT5  C (20) 


and  the  amplitude  of  the  (1  + 1 )th  oscillation  can  be  written  as: 

— °*° 
t«H  \l  1 + 


(21) 


where  Tf  is  the  period  of  oscillation. 


The  logarithm  of  the  ratio  of  the  succeeding  amplitudes  is: 


i££S 


(22) 


Illlll I I Ililiilll 


T*  = 


_ / 


?u/. 


(29) 


The  equation  of  motion  for  a small  amplitude  free-oscillation 
system  is  given  by  equation  (9)  when  the  forcing  function  vanishes, 


in 
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d*< 


f^k  NU  3 _ . f\ 

is:  t ~~  w 


(30) 


In  the  free  oscillation  technique  the  model  is  mounted  on  a flexure, 
therefore  account  must  be  taken  of  the  elastic  stiffness  in  the 
erm  and  the  flexural  damping  in  the  term. 

in  order  to  understand  the  flexural  dampint  contribution, 
an  explanation  of  the  concept  of  structural  damping3  is  in 
order.  The  equation  of  motion  could  be  rewritten  in  the 
generalized  form  as: 


+ XltA  - - UJn 


illlt. 


(31) 


The  steady  state  response  of  the  system  to  the  harmonic  excitation 
is  given  by  the  real  part  of: 


urt- 


h 


+ .2 


iurL 


(32) 


where  H(w)  is  the  complex  frequency  response.  We  can  write, 
however , 


J- 


ever,  ^ 

($$■  .2  in£f} e'* 


(33) 


where  the  phase  angle  <f>  is, 

A -4- 

Equation  (32)  can  therefore  be  rewritten  as: 

x(U>= 

which  indicates  that  the  complex  vector  describing  the  response 
x(t)  lags  behind  the  complex  vector  describing  the  excitation 
f{t)  by  the  phase  angle 


(3*0 

(35) 


XK>«iaiAlH»\c5<wL“4)=  \uiX(M 

I = COS  +■  i 318^/2.  = g1^2. 


(36) 

(37) 


Equations  (36)  and  (37)  indicate  that  the  velocity  vector  is 
90°  ahead  of  the  displacement  vector  and  W times  as  large. 
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IMAGINARY 
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FIG.  3 VECTOR  REPRESENTATION  OF  EQUATION  (38) 
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Inaepsnaenc  ox  i 
of  os cil  ia t ion . 
proportional  to 


he  race  of  strain  and  depends  on  the  amplitude 
The  energy  loss  per  cycle  was  found  to  be 
the  amplitude  squared 


(43) 


where  !■;  h 

f'  r>  A fr  »"**”*  v r*  **  h 

» a Vvi  V*  *i 


constant  of  proportionality  independent  of  the 
monic  motion. 


Similar  behavior  is  exhibited  by  any  elastic  material.  This 
type  of  damping  is  known  as  structural  damping  and  is  generally 
attributed  to  the  hysterisis  of  the  material  or  to  friction 
0$tW66  n w a o con«pO;it^  nts  of  the  system.  Even  though  this  type  of 
damping  is  independent  of  the  frequency  of  the  cycling  stress, 
it  is  a phenomenon  associated  with  cyclic  stress  as  the  energy 
dissipated  in  one  cycle  is  proportional  to  the  area  within  the 
hysterisis  loop  for  the  material  (see  Figure  3). 


Therefore,  one  may  treat  a system  subjected  to  harmonic  excitation, 
with  structural  damping,  as  a viscously  damped  system  with  an 
equivalent  viscous  camping  coefficient.  Such  a coefficient  is 
suggested  by  comparing  equations  ('12)  and  (43). 


c ^ — 

TTU) 

It  is  now  possible  to  write  the  equation  of  motion  of  a 
one-degree-of-freedcm  system  with  harmonic  excitation  and 
structural  damping  as: 

/xx  . \i/i  \ i w u<i  \ _ y*  h,  s*  I tot. 


(44) 


MKW+#5  xai  + Kxcu  = KKe1 


(45) 


Comparing  t »i“  IlOiil  ogeneous  forms  of  equations  (19)  and  (1*5) 
suggests  that: 

jf  fi/n  = X - 2 7r£lu 


(46) 


Account  must  be  taken  of  the  structural-mechanical  and  elastic 
stiffness  as  represented  by  terms  with  the  "tilda”  in  equation  (; 
Again  comparing  equation  (19)  with  psO and  equation  (30)  yields: 


- M M 

- _ ( 

Mjt 

X 

- V 

t 4 

-!*± 

+ 

X 

\ 

X 

(47) 

(48) 


FIO.  4 HYSTERESIS 


*1  ==  *•;!*  /Ui  if  r^PI«  ^ 

iii  I 


The  terns  with  the  bar  contain  the  mechanical-structural  and 
stiffness  contribution  and  the  terms  without  the  bars  contain 
the  aerodynamic  contribution. 


. 

"l TTZ 

(^9) 

— * 

IA*  ^ 

x ~ T74 

(50) 

U>*  = 

I© 


Since  we  want  the  aerodynamic  contributions  alone  we  must  evaluate 
the  mechanical  terms  by  permitting  the  model  to  oscillate  in  a 
vacuum  or  tare  condition.  Tare  values  are  indicated  by  the 
subscript  o.  The  total  damping  experienced  may  then  be 
determined  using  equations  (8),  (^6),  and  (**7)  as: 


¥■[ 


+ 2V1  ] 


(51) 


Sene  rearranging  and  combining  with  equation  (1*6)  yields: 

= irr«J  = 


which  can  be  used  to  rewrite  equation  (51). 


Solving  for  C yields: 


(52) 


WJ  / 


[^r]  = “*  WA©/uJ 4 ^ 5 

or, 

CM<t=  (5; 

where  A and  U)^  refer  to  the  values  recorded  with  the  wind 
tunnel  on  and  X*  andWj  refer  to  values  obtained  at  the  vacuus 
condition.  The  above  couid  also  be  expressed  in  terms  of  the  3d 
decrement,  g , as: 

c„,=  ^rs-<  n 


wstw  **1*^  ■ ■ ■ ■ ■ * 


*sr»f  */»  /?  *r  ; /mrs  r *?  •* 

>wBL/auu/ in  /p— D1! 


"he  aerodynamic  stiff/ies 
:he  same  procedure , 


s term  may  be  found  using  essentially 


- 1 [£fr  -’Sr]  3 -I 

*•  #•  J 


(57) 


making  use  of  the  log  decrement,  £ . as. 


id  in  coefficient  form. 


C4 


(53) 


(59) 


[*Hi  0+4r*  “ w* 


(60) 


The  equations  of  motion  of  a one-degree^of -freedom  viscously 
damped  (including  structural  damping)  system  system  have  been 
developed  from  the  general  case.  Forms  have  been  generated  which 
allow  calculation  of  the  damping  moment  coefficient,  the  pitching 
moment  coefficient  based  on  trim  angle,  and  the  slope  of  the 
pitching  moment  coefficient.  In  order  to  evaluate  these 
expressions  parameters  must  be  obtained  for  a wind  tunnel  on 
case  and  for  a vacuum  or  tare  case. 


= M sin&Jrt.  is 
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or  changing  the  variable  of  integration 


I $ a i Jt.  + M;i  $ + 

• A a i.  ’ 


-*  X 


The  steady  state  solution  of  equation  (9)  is  given  by  making 
use  of  Laplace  transforms,  equation  (9)  becomes: 


(62) 


s =.  E||fer 


(63) 


(6U) 


Now  solving  for  «C  we  get: 


ot  — 


Z 

Wn 


r^-^iis+n 

L<  <**.  1 


r^-+i 

L<  ' 


(65) 


Since  the  inverse  transform  operation^  j-f  9 is  linear,  we  can 
regard  the  inverse  transforms  separately:  The  first  term  on 

the  right  is  the  transient  solution  With  the  inverse  of  this 
tern  already  obtained  in  equation  (17).  As  has  been  pointed  out, 
equation  (18)  represents  the  angular  motion  of  the  model  during 
free  oscillation.  For  the  forced  oscillation  technique  we 
can  regard  the  system  as  starting  from  rest,  i.e.,3#*  O , 

in  which  case  the  first  term  on  the  right  c-f  equation  (655 
vanishes.  The  inverse  transform  of  the  second  term  on  the 
right  of  equation  (65)  is  as  follows: 

f s:" 

where  J 


(68) 
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Even  in  equation  (66)  we  have  a transient  solution  in  the 
second  term  which  accounts  for  bringing  the  system  from  rest 
to  the  steady-state  forcing  frequency..  This  transient  would  be 
superimposed  upon  the  transient  represented  by  equation  (17) 
if  is  not  zero.  Since  we  are  interested  in  the  steady-state 

conditions,  only  the  first  term  on  the  right  of  equation  (66) 
will  be  retained.  This  term  may  be  rewritten  in  a slightly 
different,  but  more  useful,  form  as. 


p/u>, 


2 

A 


(69) 


where 


C70) 


Resonance  occurs  when  the  amplitude  of  the  sine  function  of 
equation  (69)  attains  a maximum  value.  In  the  absence  of 
damping  resonance  occurs  when 

m-  _ , 

U>n  “ 

and  in  the  presence  of  damping  resonance  occurs  when 

uM-  1 ^ 

U >a 


- J\~* 


(71) 


(72) 


Resonance  is  often  said  to  occur  when  the  forcing  frequency,  Wf*, 
equals  the  undamped  natural  frequency,^.  This  statement  is  true 
only  when  the  system  has  no  damping,  although  it  may  be  an 
acceptable  approximation  for  low  levels  of  damping. 

The  above  expression  is  easily  derived  by  taking  the  derivative 
of  the  amplitude  in  equation  (69)  with  respect  to  (wf/u>n)  , 
equating  this  derivative  to  zero  and  solving  to  obtain  equation 
(7?).  Inserting  (72)  into  (69)  gives  the  amplitude  at 
resonance  as 


KC^ 


1 \ ?/u£  _ v 

xj v-fy*  = 


(73) 
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Also  at  resonance,  the  phase  angle,  (,eq.  (70))  becomes. 


(7*0 


The  angle, i p,  is  known  as  the  phase  angle  and  it  may  be  seen 
from  equation  (69)  that  the  angular  displacement  lags  the 
applied  moment  by  the  angle  ip.  The  applied  moment  from 
equation  (*0  is: 


fACtV=  fASwWjt 


(75) 


and  the  angular  displacement  follows  from  equation  (69)  as, 

<*  (t^  = K Sin  - ^ ) 


(76) 


The  relative  position  of  the  moment  and  angular  displacement 
functions,  M(t)  ando<(t),  respectively,  are  given  in  Figure  5. 

Since  we  are  considering  a lightly  damped  system  (damping 
ratio  less  than  0.1)  a first  guess  at  the  phase  angle  might  be 
90  degrees.  However,  the  damping  ratio,  though  small,  is 
influential  at  resonance  both  in  the  angular  displacement  amplitude 
CviZf,  eq.  (73))  and  in  setting  the  phase  angle  exactly. 

Clearly  equation  ( 7 *0  may  be  rewritten  as* 


swop  = = / - T^r 


(77) 


Since  the  actual  phase  angle  will  be  slightly  less  than  ir/2, 
we  may  introduce  a small  positive  quantity,  e,  as 


Equation  (77)  becomes 


si* ofr-O3  Cos'e=  I-  Tt* 


(78) 


(79) 


or 


Si<\£^ 


0-f) 


*\Y*. 


or  solving  the  above  expression  for  the  damping  ratio,  C,  we  get 

_ Sif\  £ 


[UW'c'l 


vz 


(80) 


2H 


FIG.  S MOMENT  AND  ANGULAR  DISPLACEMENT  VS  TIME 
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Clearly,  equation  (79)  allows  the  following  approximation  for 
small  values  of  damping  ratio. 


£ = f 


(81) 


Thus  in  order  to  measure  the  damping  ratio  (and  hence  the 
damping-in-piteh  derivative  from  equation  (15)) it  is  clear 
that  the  phase  angle,  !/>,  or  the  term,  e,  must  be  measured 
accurately.  Equation  (69)  together  with  the  first  and  second 
derivative  may  be  written  as. 


iz  - w>4.K  Cos 


(82) 

(83) 

W) 


Inserting  the  above  expression  into  equation  (62)  where 
appropriate  gives, 

- X UJj.  K 5 C-os(«Ujt-V)Sin(u^t-y)  Jt  - ^tos(Ht-*)Jt 

- )Cos(w^-f)Jt  = (85) 

S^Sin  u»ftCo&C«Ht- ^)dt. 

Clearly  the  first  and  third  terms  on  the  left  are  zero  so  we 
are  left  with, 

.'T 


• y 

- ft*  0>4  K*  f)<tt  Sin  «^t  CoiWjfc 

Cos*  it  + MMfcfen 

o 

The  first  term  on  the  right  is  easily  shown  to  be  zero  so 
equation  (86)  becomes, 

MKu^SinY  £ 


(86) 


or 


(87) 
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It  was  pointed  out  in  the  case  of  free  oscillation  (equations 
(30)  and  (3D)  that  structural  damping  must  be  removed.  Thus, 
equation  (87)  becomes. 


or 


h/LSinVo 

V1* 


) 


(88) 


The  above  expression  is  considered  to  be  the  basic  data-reduction 
equation  for  measuring  the  damping-in-pitch  derivative  using 
the  forced  oscillation  technique.  This  equation  appears  as 
equation  (6)  in  reference  8,  as  equation  (42)  in  reference  2 
and  as  equation  (25)  in  reference  9.8*9  Along  the  lines  of 
equation  (56),  equation  (88)  might  be  written  in  coefficient 
form  as 


c = - gv.rtASi«V_  M0S.nV5,l 

M1  9S41  L H K <uf  K.  i 


;89) 


It  should  be  pointed  out  that  for  a finite  amount  of  damping 
resonance  is  reached  somewhat  before  ajf/a)n  equals  unity  (.see 
equation  (72))  and  at  resonance  the  phase  angle  is  less  than 
90  degrees  (see  equation  (74)).  From  equations  (72)  and  (74) 
we  have  for  effects  of  the  order 


“H/ib,  = i- 

SrtvT|»R=  i-Ys-£* 


(90) 


(91) 


8 

Mackapetris,  L.  J.,  "A  Forced  Oscillation  System  for 
Measuring  Damping  Derivatives  at  Subsonic  and  Transonic  Speeds,” 
Naval  Ship  Research  and  Development  Center  Report  2627, 

November  1967 

q 

Wiley,  H.  G.,  "A  Method  for  Accurately  Measuring  Dynamic 
Stability  Derivatives  in  Transonic  and  Supersonic  Wind 
Tunnels,"  Presented  to  AGARD  Specialists  Meeting,  Aircraft 
Stability  and  Control  (Available  from  NASA,  Langley  Research 
Center),  April  1961 


27 


NSWC/WOL/TR  75-84 


Quite  obviously  for  lightly  damped  systems  (c  about  0.1  or  less) 
there  is  negligible  error  in  oscillating  the  model  at  the  undamped 
natural  frequency. 

Actually,  the  undamped  natural  frequency,  u^,  is  a somewhat 
artificial  quantity  in  that  it  is  not  directly  raeasureable 
damped  natural  frequency,  a)d>  and  the  damping  ratio,  C, 
which  is  expressible  in  terms  of  the  log-decrement,  6 
(equations  (22 )— (29) ) . From  equation  (17): 


*\Sz 


(92) 


From  equation  (72),  the  ratit  of  the  forcing  frequency,  wf,  to 
the  undamped  natural  frequency,  wn,  under  conditions  of  resonance 


is. 


4\ViL 


which  gives. 


A* 

«*._  JJzlLL 

*sv  ci-pr- 


(93) 


(9*) 


where  the  relationship  with  the  phase  angle  follows  from  equation 
(74).  Thus,  at  resonance  the  ratio  of  the  driving  frequency, 
uf,  to  the  damped  natural  frequency,  u>d,  is  equal  to  sint{»R. 

While  the  primary  goal  of  the  small  amplitude  forced  oscillation 
technique  is  to  obtain  the  damping-in-pitch  derivative, 

it  is  possible  to  also  obtain  the  local  tangent  to  the  static 
pitching  moment,  i.e.,  From  equation  (69)  we  have. 


and  from  equation  (70) 

CoSTj»  = 


[o-(^)Vf(£ir 

i-  fe') 


(95) 


(96) 
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m 1 


Thus  equation  (95)  may  be  rewritten  as. 


V/Sifc  -f  \ - P CosTp 

< (I- 


(97) 


Nov  from  equation  (48) 


1/ /Sit  - P^OS*#  _ pCoS"^ 

ft\uun  ]"  — ; T — TT . "af  z. 

* u>*-u>£  -BL,^ 


or 


- X ' “x" 


- W?*  s 


(98) 


(99) 


In  a vacuum  the  structural  damping  contributions  are  alone  present. 


M = - McCoslft,  _ j w 


K. 


(ICO) 


Using  equation  (100)  In  equation  (99),  the  static  moment 
derivative  due  to  aerodynamic  effects  is 

- r»A _ MpCoS  (u^- <o£ ^ 


aoi) 


The  assumption  has  been  made  that  the  mechanical  damping  is 
proportional  to  the  oscillatory  frequency.  Vacuum  conditions  are 
designated  by  the  subscript  "o."  Equation  (101)  may  be 
rewritten  in  coefficient  from  as 


C = - 


Hi 


\_ 

9*4 


(102) 


Equations  (101)  and  (102)  are  the  data-redueticn  equations  for 
the  determination  of  the  static  moment  from  the  small  amplitude 
forced  oscillation  technique.  Equations  (101)  and  (102)  appear 
as  equation  (24)  in  reference  8,  as  equation  (43)  in  reference  2, 
and  equation  (26)  in  reference  9* 

Further  simplifications  are  possible  in  equations  (101)  and 
(102).  If  the  drive  frequency  is  identical  in  both  the  wind-on 
and  wind-off  conditions  then  the  second  term  in  equations  (101) 
and  (102)  vanishes  to  give. 


r _ _ J ( >ACes-»_  _ M„Cos 

<*sA  ( K K«  ) 


(103) 
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In  driving  at  resonance  the  phase  angle, 

Co STPe  - - « t 


is  such  that 


(10ft) 


a result  that  follows  readily  from  equation  (7ft).  Clearly 
equation  (10ft)  shows  that  for  the  lightly  damped  system  cos'f'R 
will  be  very  small  and  approximately  equal  to  the  damping  ratio. 


Another  concept  worth  noting  is  the  ”Qrt  or  quality  factor  of 
the  system.  The  Q of  the  system  is  defined  as  the  ratio  of  the 
amplitude  under  dynamic  conditions  to  the  amplitude  under  static 
conditions,  i.e., 

iaio 

as  (us) 


The  static  deflection, , may  be  obtained  from  equation  (19) 
by  setting < = «£  = 0 to  obtain. 


(106) 


The  appearance  of  frequency, in  the  evaluation  of  a static 
effect  might  be  misleading.  In  equation  Cl6)  it  is  shown  that 
the  undamped  natural  frequency,  wn,  is  determined  by  the  static 
moment  derivative,  From  equations  (76)  and  (95)  we  get  for  Q, 


Q 


|3t*ai  _ 


(107) 


For  osci 


resonance  we  may 


find  Qp  using  equation  (72)  as 


J -v  J— 

rfQ-V)  ~ *7 


(108) 


or  using  equation  (23) 

<?« 


3L 

6 


(109) 


Equation  (107)  indicates  that  for  a lightly  damped  system  it 
is  necessary  to  oscillate  almost  exactly  at  resonance  to  see  any 
amplification  of  the  static  response.  Equations  (108)  and  (109) 
show  that  when  resonance  is  closely  approached  there  will  be  a 
sudden  increase  in  the  amplitude  of  the  model.  Equation  (70) 
points  out  that  for  a lightly  damped  system,  the  phase  angle  will 
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be  zero  until  very  close  to  resonance  at  which  point  it  will 
be  near  to  ti/2  (although  as  equations  (73)  and  (8l)  indicate  a 
more  precise  value  would  be  ti/2  - 5).  The  consequence  of  making 
dynamic  measurements  on  a high  Q system  is  that  there  might  be 
some  difficulty  in  practice  in  ’treasuring  the  term  Msinift 
accurately. 


FREE  OSCILLATION  DATA  REDUCTION  PROCEDURE 

As  an  introduction  to  the  data  reduction  technique  a discussion 
of  the  principle  of  least  squares  (Ref.  *i)  and  the  Method  of 
Differential  Corrections  (Ref.  5)  is  in  order. 


The  principle  of  least  squares  in  simple  terms  states  that 
the  curve  which  best  represents  a set  of  data  points  is  the  one 
for  which  the  sum  of  the  squares  of  the  residuals  has  a minimum 
value.  The  term  "residual”  as  referred  to  in  the  previous 
sentence  implies  a difference  between,  for  a given  obscissa,  a 
data  ordinate  and  the  curve  ordinate  (see  Figure  6). 

If  a function  y = f(x,a,b,.  . .)  is  to  be  fitted  to  a set 
of  data  points  (Xj,  Y3J,  the  major  problem  involves  finding  values 
of  the  constant  terms  a,  b,.  . . which  will  generate  the  function 
which  will  best  represent  the  given  data.  Application  of  the 
least  squares  principle  requires  determining  the  values  of 
a,b,.  . . which  will  make  the  sum  of  the  residuals  squared, 

*V?  , be  a minimum.  This  may  be  accomplished  by  developing 

the  set  of  simultaneously  solvable  eauations. 


4 <x 
* *> 


= O 


- o 


(110) 


In  the  function  to  be  fit  f(x,a,b,.  . . ) is  linear  in  a,b,.  . ., 
then  the  residual  equation  (110)  will  also  be  linear  and  the 
system  may  be  solved  with  little  difficulty.  In  the  functions 
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which  must  oe  used  to  represent  the  wind  tunnel  motions 
presently  being  considered,  the  parameters  a,b,.  . . are 
arguments  of  exponentials,  sines,  and  cosines.  Under  this 
requirement  the  residual  equations  become  non-linear  and 
cannot  be  directly  solved.  The  Method  of  Differential 
Corrections , which  will  be  discussed  next,  provides  an  iterative 
method  of  solution  to  any  degree  of  accuracy. 

Given  a set  of  N data  points  (*'*i ,Yi ) and  a function  of 
y± =f  (:<±  ,a,b  ,c)  which  is  to  be  fitted  to  the  data. 


The  residuals  are  given  by. 


Vj  =.  - Y; 


(in) 


Let  the  parameters  a,  b,  and  c be  defined  as  follows, 

ol - a0  * 
b = bo  4-  & Y> 

Cl  = C0  + A C. 


(112) 


where  aQ,  bD,  and  c0  are  first  approximations  of  a,  b,  and  c. 
The  problem  is^not  to  determine  the  corrections  £a,  £b,  and 
Lc  so  that  Zv-s  2 is  a minimum. 


Equation  (111)  may  now  be  written. 


(113) 


Expanding  by  Taylor’s  theorem. 


Yi+Vi  -4  (x*Ma,,bo,C^4- A*C**V*<0  +• 

AV>  + ^P^Co’)  4 OrJtr 

where  f.  means  f(x. ,a  ,b  ,c  ). 

X jl.  o o c 

If  the  first  approximations  are  reasonably  good,  the  &*s 
•will  be  small  and  the  higher  order  terms  may  be  dropped. 


Let  us  define  another  residual 


♦•he  difference  between 


the  ordinates  of  the  first  approximation  curve  and  the 
data  points, 

ft;  = X;  - 4 Cx-Ma.,Y>.,  c.') 


Equation  (11*0  may  now  be  written. 


+*C(&V 


(lib) 
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i-’U,-  - It  / *>—  b U 


The 

condition  for  £y . ^ 

?A  hp  ^ p 

w v[w  a «• 

.»o,**'=o 

a&a 

J j b c. 

ubstituting  (115)  i 


U. , B; , , C;  = ° ’ 

= ys-a^(^o-A;+M>»B;4.Ac.Ci-R,-i'=  O 

= Z£ A\  + O 

-jTXST  = C^«.Ai+W»6-.thcCi-R;^  =o 

51  * *c,  ( **  A.+ m.%.+ t>tc, -*0  - o 


(117) 


In  matrix  form, 

r^.  7. 


£A.  ^A.B;£h;C. 


^B«A«  £B>  £ B.c. 

^C;Ai  £C;B|.  £C.x. 


fcV> 

e»c_ 


£A;R? 

^BiR; 

*C;R. 


'iting 


(113) 


Equation  (113)  are  known  as  the  normal  Eouatio^s*  their 

luticn  results  th<=-  vai»e-  * 7 ’ 

— • — -*•-  /aiues  Oi  na,  iio,  ana  «e. 


atx ting  eauation 
equation  (118)  as  a g 
right  by  the  column  c 


;r*ix  c 

SOSO  j 


a matrix  02  coefficients  and 
coefficients  augmented  on  the 

A t"  ■? 

~ * *“■  9 


c,  *J 

Cx  Rx 


(119) 


£8,*  A. 

£C;A; 


£A.B,  £A .«? 

£B;c;  * ftR, 
£C;%;  £C>  £ C;R. 


(120) 
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A simple  algorithm  by  means  of  which  matrix  (120)  may  be 
formed  directly  from  matrix  (119)  is  evident:  To  form  the  first 

rev;  of  (120)  multiply  each  row  of  (119)  by  its  own  first  term 
and  sun  the  columns;  to  form  the  second  row  of  (120),  multiply 
each  row  of  (119)  by  its  own  second  term  and  sum  the  columns; 
to  form  the  third  row  of  (120),  multiply  each  row  of  (119)  by 
its  own  third  term  and  sum  the  columns. 

The  residual  equation  matrix  (119)  could  be  very  large,  since 
there  is  a row  for  each  data  point.  Examination  of  t.*e  above 
algorithm  shows  that  it  is  not  necessary  to  generate  (119) 
as  a matrix.  Ke  may  generate  the  first  row  of  (119),  form  the 
necessary  products,  and  add  them  to  the  terms  of  (120) 
immediately.  The  first  row  of  (119)  may  then  be  replaced  by 
the  second  row,  the  necessary  products  formed,  ana  added  to 
(120);  repetition  of  this  process  N times  produces  the  complete 
normal  equation  matrix  (12*0). 

The  solution  of  (120)  provides  values  of  Aa,  Afc,  and  Ac. 

These  are  used  to  correct  the  approximate  parameters  aQ,  bQ, 
and  cc;  these  corrected  parameters  are  then  taken  as  improved 
approximations  and  the  whole  procedure  repeated.  The  iteration 
is  continued  until  there  ceases  to  be  significant  change  in  the 
parameters;  the  process  is  then  said  to  have  converged. 

It  is  more  convenient  to  test  for  convergence  on  a single 
quantity  than  on  the  several  parameters.  IRj2  Is  a suitable 
quantity;  -when  this  ceases  to  change  significantly,  the  fit 
is  as  sc-od  as  can  be  achieved. 


iquaticn  (13)  may  be  rewritten  in  the  general  form, 

oC  - K C + (°*t 


(121) 


where  the  sine  and  cosine  functions  differ  by  a phase  angle  ana 


5 f*  =*^11  r= 


^ scGulil Zion  svst^oni 


the  term  a~+cr;  account  for  trim  and  misalignment  of  the  absissa 
of  the  data  acquisition  system  and  the  velocity  vector.  In 
the  fitting  procedure  arr-.+oy  cannot  be  determined  separately 
so  a single  factor  K-3=  (d-p+ay  ) Is  determined.  Equation  (121) 
is  the  fitting  model' which  is  utilised  to  represent  the  one- 


The  independent  variable  of  the  fitting  model  is  time. 
The  data  is  assumed  to  consist  of  corresponding  values  of 
missile  attitude  angle  and  time.  The  method  of  differential 
corrections  imposes  no  Inherent  requrement  for  a constant 
time  Interval.  Data,  however,  is  usually  presented  In  terms 
of  a constant  time  Interval,  and  great  simplification  of 


logic  results  whs 
uniform  scale. 


pc  t ir  r t ns 


lone  with  resoect  to  a 
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K and  g are  the  magnitude  and  orientation  of  a vector  ana 
at  t=0.  If  a section  of  data  late  in  a flight  or  wind  tunnel 
run  is  fitted  in  terras  of  real  time,  K and  g will  he  extrapolated 
back  to  zero  time  in  accordance  with  the  current  values  of  A and 
u.  It  is  therefore  incident  that  the  time  array  should  be 
normalized  so  that  zero  time  is  located  within  the  section  of 
data  being  fitted. 

wind  tunnel  data  are  always  at  least  slightly  non-linear; 
that  is,  a and  w change  somewhat  with  time.  The  values  of 
A and  a obtained  by  fitting  a linear  model  to  such  data  are  mean 
values  related  to  the  mean  or  middle  time  of  the  section  of 
data  being  fitted.  Since  K and  g are  related  to  zero  time,  the 
time  array  should  be  normalized  so  that  zero  time  is  at  the 
middle  of  the  section  of  data  being  fitted;  all  parameters  will 
then  be  obtained  related  to  the  same  instant  of  time. 

The  parameters  to  be  determined  by  fitting  the  model  to  the 
data  are  K,  A,  w,  $,  and  K3.  Reasonably  good  first  approximations 
of  these  quantities  are  necessary  to  start  the  differential 
corrections  procedure.  The  final  results  do  not  depend  upon  the 
accuracy  of  these  approximations;  the  only  requirement  is  that 
they  be  good  enough  to  result  In  a convergent  iteration.  These 
first  approximations  are  illustrated  in  Figure  7. 

K-2  may  be  determined  as  the  mean  of  the  two  extreme  points  of 
minimum  amplitude. 

fn— ^ '1 

u may  be  determined  from,  u = ------  3,  where  n is  the  number  of 

extreme  points  and  tn  Is  the  time  interval  between  the  first  and 
last  of  the  extreme  points* 

g may  be  determined  as  g * wt$,  where  Is  the  time  interval 

between  the  normalized  time  zero  (middle  point  of  section  of  data 
being  fitted)  and  the  preceding  positive  maximum.  A poor  first 
approximation  may  result  In  a negative  value  for  K;  the  absolute 
value  will  be  correct,  however. 

K Is  determined  as  the  distance  from  the  Kq  line  to  the 
Intercept  at  normalized  time  zero  of  the  envelope  of  positive 
maximum  points. 

A*0  will  usually  provide  a reasonable  first  guess  for 
the  A term. 


The  maximum  allowable  error  in  the  approximation  of  & depends 
upon  the  resulting  maximum  phase  difference  between  the  model  and 
the  data.  For  a given  error  in  u,  the  phase  difference  Is 
proportional  to  the  number  of  cycles  included  in  the  section  of 
data  being  fitted;  therefore,  greater  accuracy  is  required  when 
fitting  a long  section  of  data. 


G.  7 FIRST  APPROXIMATIONS 
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The  derivatives  required  to  form  the  residual  equations  can 
now  be  expressed  as. 


6;  - _ -.It.' 

" ' v 


CoS(Wfc;+<$) 


Ci  = Kt;  B; 

Di  = ^ - Ke^smtu*.;**) 

Ei  “ trcr  - 

R;  « •*;  - 

The  normal  equation  matrix  differs  from  equation  (120)  only 
in  the  addition  of  two  rows  and  columns  for  the  additional 
derivatives  D and  E.  The  summations  of  equation  (120)  are 
for  i=l  to  N where  N is  the  number  of  data  points  to  be  fitted. 


Following  is  a listing  of  a FORTRAN  IV  subroutine  which  has 
been  used  to  fit  the  data  of  the  current  discussion. 


* ■ > - > PLUMBUM 


NSWC/WOL/TR  75-84 


FPOGfiAP  CAHDAKP(INPUT«£lTPL1  » 1 APE  5 = 1 **HJ  I * 1 AF  Et.=CU  I Pl,T  • T API-  /t  • 
1TaPE1*TAPE?*7APE9S) 

DIMENSION  CO(Fl*CC»l>  *1  «LN(*>C)  *TI*h  (Midi)  .«  <F*-*i) 

PE  AH  <5*  10)  KM  tx 

IE  (KM.'K  ,\f  .0)sh  At.  (5*  J (• ) <LM,f.  | ] ) * 1 = i .hn;|  ) 

Id  FOPEAT  (?<«n) 

CALI.  FU1  I(CO*?OCO»»x<*(M'«r  Hi  K.I.i-uM 

Oc  500  «*7  = 1 • 2 

PEAL' (5*11)  TrFG.TfcM) 

11  HWAT  C2F1P.6) 

mF«D  ( 1 1 IMM3*NS»lSA«IfcA«fxA*j|-r.>S*jFHS2*Jh<LSZ*txMjK 

IF  (JFPv.S.OT.ftOOO)  ST  CPI 

PE  AH  Cl) ( (C0(L*K) *K=l«r  A) *L  = 1,JFPPS) 

jFkF  S = JF  FPS-1 

LCW=0 

li-l»=l 

i;c  55  LM.JhF*'? 

IF(L.EO.l)  TI*F<L>=0.0 
IF(L.C-T.l)  1 1 >-‘E  (I. ) = I 1 f*F  (L-l ) ♦• 004 

IF  ( (T  IN*F  (L  ) .(  T.T-tP)  »«  i.  (T  (L)  • fcl.IrMi)  ) C-C  TC  55 
LCN=LCN*1 

IFdCK.EC.l)  TIM;  U CM=<).0 

IF  (LCk.GT.1 ) TIM (LC^)sTIFE (LI  \-l)*.G04 

Qf'  (LCK  • l ) =CO  (L  * 1 ) 

CS  CONTINUE 
JF  KPS=LCK 

CALL  VALLE  (NPI.N*  JFP^5*T  IKf  *CO*  a • I * f.  *CM  « f f-F  0 • 1 1 M>  > 

JFMS  = <*00 

CAI_L  FIT  (NkUK*  JFPKS«T  IM.  *CC  *1  L*  • a *CM  ) 

SOU  CONTINUE 
STOP 
h NO 
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SUBROUTINE  VALUE  <NWUN..,jFHFS.TlRt  „C C • a ♦ ! *U* CM  ♦ Tc*EG . TEMJ ) 

DIMENSION  CO(f  OCO«1 ) *1 1 f- £ (MOO*  >=kAX  (?M  • I FAX  ittb)  • P»lM(?b)  vIPlMirE 
1 ) •A(2f-»3) 

LKlN=l 

IRAX=1 
X*IN=0.0 
XB AX-0 • 0 

1)0  10  w = ?«JFf*/S 

IF(CC(J«l).n*CC<  *.-!♦!>)  = 1 ,<• 

IF  (X* AX. EC. 1.0)  F^AXILF/ *)=C0(J-1 «1) 

IF  <Xi* AX. EC.  1.0)  TF'AXtL^AXJsl  lf'h(c-l) 

IF (XNAX.EC.1.0)  LNAX=1*AX>1 

IF  (LRAX.EG  • 1 2 ) GO  TO  IS 

IF(J.F0.2>  LKAX=I 

IF  (X* AX. EC.  1.0)  aMN=O.0 

IF  (CC(J*1)  .GT.CC  L-1»U)  X*Ik=A«IN*1.0 

IF  (X* IN. EC. 1.0)  Ff'IML»'IM=CO(J-l.l) 

IF  (XMIN.EC.l .0)  TMN Cl.^lN)  =T  IKE  <»;-l) 

IF  <X* I*. EC. 1.0)  LKIK=I.MNM 
lF(J.F<i.2)  LR  IN’=  1 
IF  IX"IN.FC.]  ..»)  XFA**O.C 
in  CONTINUE 

IE.  IMIxC.tC.n  .-mTfc(f.FS)  l<*LN 
c 9 FORMAT  (lhl.lOX*<»HMJN  *I?.10X.M^*-INC  UEF//) 

IF  <I»n.EC.2)  w^I  TE  N«LN 

c»  FONFATClUf  10**4hH.«v  •IS*lO>»7n^IT*C  ON//) 

► *»  1 T 2 (6«1C1 ) THfcf-.IEM 

teX'  I T E (F« i CO)  ( cm/,  f ( I)  *1  = 1 *M 

WHITE  (F..  101)  UMAX  (I)  *T  = I.‘*) 

•FITE  If. ICO)  (PKIN(I) *I  = lfS) 

VFITF  (F.,1C1 ) (TMIMI)  *1-1*^) 

100  FTKVAT ( 1 OF  1 0 • 1 ) 

101  FORM!  (10F10.3) 

CNT=TMAX(1) 

A<***1  JsCPPA*  <1  ) «R*IM1 ) )/2.f> 

A (A  • 1 ) =1  •?  71 

FPE(.s(TRAX<  9)*1)P«(l))/(j,ft 
A ( 3* 1 ) =6 .2P31 P/F  REQ 
A(?.1)=FPAX(1)-A<S»1) 

*M.|M  = PMAX  ( U -l-KAX  (V) 

UF.NP  = (TP AX  ( <•)  - TWAX  (1  ) ) « (F'-AX  <£)-*  (h«l  ) ) 

ACl*l)=AMM/rFfjfc 

WHITE  ((.,50  E.CUN*  (A(M)«I  = NN)«(.M 
*:0  FORMAT  (IS«6F10.3) 

RETURN 

END 
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SUBROUTINE  F IT  (K  Wt**  JF  *FS  • T I M *CC*LCN*A,CNT) 

DIMENSION  NC5)  *X  (<)*<.0OC)  » A (26*3 I *CC  (<*(!00*1)  *T  IKE  UOOO) 

EXTERNAL  VF 
NPAR=«* 

N (1 )=NPAR 
N(2)=jFF*S 
N(3>=2 
N(*)=l 
K (5) =0 

DO  10  1=1 ,NPAfi 

10  A(I,3)=I 

DC  20  L=1 , JFFMS 

X <1,L)=TIRE  (L  ) “CM 
X(2*L)=CC(L*1) 

20  CONTINUE 

DFLT=1.0F-6 

DC  777  NPASS=1*S0 
1F{KFASS.EC.1?>  RETURN 
CALL  LSCSLU (N  * X • A • VF  *OFLT ) 

IFCA<?**1).G1,0.)  GO  TO  62 

*HTE(6*60)  (MI*1)*I  = 1*KFAR)«A(?G*1)«A(2£*2)*a(26*3) 

60  FCHKATt  *F 1A.7/5F1A.7/) 

IF  ( A C 1 « 1 » .LT.O.C)  rfturn 

777  CONTINUE 

f?  *PITE(fc.60>  (A(I«1)«1=1 *KFAh ) * A ( 26  • 1 ) • A (2b  *2 ) * A c 26 * ? ) * 

I (A(l*2>  *I  = 1*NFAF> 

K (2)=NlNO(LCN**COO> 

JFWRS=HNO(LCN*<iUOO) 

DC.  11  1 = 1 .NPAR 

11  A « I *3) -1 

DO  21  L = 1 » JFPV.S 
X(1*L)=T1RF (L)-CNT 
X(2*L)=CC(L*1> 

21  CONTINUE  ’ 

DC  778  NFASS=1 *50 
IF  (NFASS.EQ.12)  RETURN 
CALL  LSCSUB(N*X*A*VF tDELTl 
IF (A(?6*1).GT.0.>  GO  TC  63 

tePlTE (6*60)  (A(I*1)«I  = 1 *UPAk) * A (26* 1 ) * A (26*2 ) *A(2t>*3) 

IF  (A(l*l)  .LT.O.C)  RETURN 

778  CCNTTKiE 

63  %*F«  I T E (6*60  (A  (1  .1)  *1  = 1 *NPAF ) • A (2t  * J ) * A (26*2)  *A  (26*3)  * 

1 ( A ( I •?) ,1=1 *NF Ah ) 

65  DC  25  K=1 , JFPKS 

C0(F*1)=EXP(-A(1*1)*(TIRE  (K)-CM  I )*A(2,1)*SIK(A(3*|)*(IIMF  (K)-CM) 
1 ♦ A (A  « )’)  ) *A  (5) 
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25  CONTINUE 
SL'Ml  = 0.0 
K^ITE (6*200) 

200  FORMAT (//•  RLN  IH(  TInE-CM.  INPl 

1 L I F F * ) 

DO  30  K=1 * JFPm$ 

DIFF=CO(K.l  )-X  (?.M 
SL«l=SUFl*DIFF*«2 

30  CONTINUE 

D£V=SCKT (SUM1/JFRKS) 
teKITE  (6,101)  CEV 
101  FCRMT (E14.7) 

DC  31  K=1.JFRMS*20 
DIFF=CO(K,l)-X(2»K) 

TIK=TJME (tO-CNT 

teHTE  (6*100)  NhLN*TI('E(N)  *TIf-*X(2,K)  *CO(K*l)  ,DIFF 
100  FORMAT (15, ?F 12*3, 3F 14*7) 

31  CCNTIFUE 
RETURN 

END 

END 


j»»M**i«rtw«iMiiii<iiiiii  1 1 liiMi  mrrniiiTiiiHIttilfi'n'iir1  11 " ' l****MI*>11” 
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SUBROUTINE  VF(N«X*DF«A) 

DIMENSION  OF  (i)«X(l)«A(l)*N(l) 

EOLATION  TO  FIT.***  F=t  XP  (-A  ( ] ) «xl«A  (2)  *SI  tv(A  (•}  *A*A  (4)  ) *A  (5) 

= TI*E-SHH  CCr%ST*i«T 

SA3=SIN'(A  C3)«X(])«A(4! > 

EAl=EXP(-A(J)«xm  » 

OOH  =EA1«A(2)*SA3 
OF  I2)=EA1«SA3 
OF  ( 1 ) =-X  Cl) «DLM 

SIM  =A(2)*EA1*C0S(A(3)«X(1)«A(4)) 

OF (3) =X (1)*SUP 
OF (4>=SIP 
OF  (51  = 1. 

X (9) =DUP*A (5) 

RETURN 

END 
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THE  SMALL- AMPLITUDE  FREE-OSC ILLATIOIJ  SYSTEM 


The  Supersonic  Tunnel  Association  standard  ten  degree  cone  w as 
selected  as  the  test  configuration.  The  model  has  interchangable 
nose  sections  which  resulted  in  bluntnesses  of  .0167,  .1,  and  .25 
when  moted  with  the  conical  afterbody  which  had  a 4.5  inch  diameter 
base.  Several  ballast  rings  are  located  in  the  model  near  its  base 
to  ensure  proper  static  balancing.  The  forward  end  of  the  cone  was 
fit  with  a ring  holder  which  attached  the  cone  to  the  flexure  mount 
of  the  supporting  sting.  The  cone,  noses,  flexure  ring  holder,  and 
a sample  ballast  ring  are  shown  in  Figure  8. 

The  basic  components  of  the  sting  support  and  model  system  are 
illustrated  in  Figure  9.  The  sting  is  designed  to  permit  the 
attachment  of  seven  comraerically  available  torsional  flexure  mounts 
to  allow  for  a variety  of  torsional  stiffness  and  load  carrying 
characteristics,  which  may  be  required  to  meet  a specific  test 
objective.  For  the  present  test  5/8  inch  cross  flexures  of  spring 
constants  of  K=106.0  and  13.3  inch-lbs/radian  were  selected. 

Using  flexures  of  differing  spring  constants  allows  evaluation 
of  frequency  effects. 

The  model-sting  system  is  injected  Into  the  flow  field  at  the 
desired  sting  angle-of-attack.  Secondary  oscillation  of  the 
support  system  initially  appeared  to  be  a problem.  To  alleviate 
this  problem  the  angle-of-attack  servo  mechanism  is  clamped  to  the 
wind  tunnel  superstructure.  An  initial  angle-of-attack  and  the 
subsequent  resulting  small  oscillations  about  the  stirg  angle-of- 
attack  is  achieved  by  a pneumatically  actuated  piston,  rod,  and  tripping 
lever  system.  Application  of  air  pressure  to  the  rear  of  the 
piston  drives  the  piston  and  rod  forward  causing  the  tripping 
lever  to  rotate  outward  and  contact  the  model.  The  model  is  caused 
to  rotate  about  the  flexural  pivot  point,  with  the  maximum  deflection 
angle  being  dependent  on  the  height  of  the  lever  hammer  head. 

As  the  piston  moves  further  forward,  the  level  Is  quickly  released 
which  allows  the  model  to  freely  oscillate  about  the  pivot  point. 
Application  of  air  pressure  on  the  other  side  of  the  piston 
retracts  the  actuating  rod  and  sets  the  system  for  another 
triggering  sequence.  Several  trips  may  be  achieved  during  each 
wind  tunnel  run,  allowing  several  damping  records  to  be  obtained. 

A data  record  is  also  obtained  prior  to  and  just  after  a wind 
tunnel  run  to  allow  determination  of  tare  damping  values  under 
near  vacuum  conditions. 


system 

n 

tigated 

thermal 


gages  on  the  flexural  members.  In  the  present  application,  the 
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FIG.  10  SCHEMATIC  OF  SIGNAL  PATH 


PIG.  11-  Model  and  Support  Sting  Mounted  in  Hypervelocity  Research  Tunnel 
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FORCE  OSCILLATION  SYSTEM 

The  forced-oscillatl on  technique  is  based  on  the  principle 
that  when  a simple  mechanical  system  is  externally  forced  to 
oscillate  near  its  undamped  natural  frequency,*  the  input  force 
(or  moment)  required  to  maintain  the  motion  is  just  that  necessary 
to  overcome  the  inherent  damping  in  the  system  (Ref.  2).  Hence 
a measurement  of  the  forcing  function  is  equivalent  to  a measure- 
ment of  the  system  damping.  In  theory,  the  oscillation  need  not 
be  exactly  at  the  resonant  frequency,  as  a measurement  of  the 
forcing  function  and  phase  difference  between  the  input  force 
(or  moment)  and  the  system  response  is  sufficient  to  provide  a 
means  for  deducing  the  damping.  However,  in  practice,  the 
technique  becomes  less  useful  as  the  oscillation  frequency 
departs  from  the  system  natural  frequency,  and  most  forced- 
oscillation  damping  measurement  systems  operate  at  or  near 
resonant  conditions. 

The  forced-oscillation  mechanism  used  is  illustrated  in 
simplified  from  in  Figure  12.  The  model  is  attached  to  the  pivoting 
head,  which  is  forced  to  perform  a single-degree-of- freedom 
oscillation  about  the  pivot  point  by  a yoke  and  eccentric  crank 
attached  to  a rotating  shaft  driven  by  a one-third  horsepower 
universal  motor  located  in  the  drive-unit  housing,  as  indicated. 

The  balance,  designed  and  built  at  the  NASA  Langley  Research 
Center  for  a hypersonic  free-jet  wind  tunnel,  employs  the  90-degree 
drive  shaft  turn  in  order  to  minimize  undesirable  tunnel  blockage 
and  heating  effects  which  would  result  from  a motor 
housing  mount  located  on  the  sting  axis.  Oscillation  frequencies 
from  about  one  to  50  hertz  can  be  obtained  and  oscillation  ampli- 
tudes up  to  two  degrees  can  be  generated  by  the  currently 
available  components.  Operation  of  the  balance  is  restricted  to 
maintaining  the  plane  of  oscillation  coincident  with  a vertical 
plane,  although  the  balance  may  be  pitched  either  up  or  down  to 
a maximum  angle  of  30  degrees  and  may  be  yawed  to  angles  restricted 
only  by  load  limitations  on  the  balance  head.  Hence,  in-plane 
damping  is  achieved  by  pitching  the  balance  in  the  oscillation 
plane,  and  out-of-plane  damping  is  obtained  by  yawing  the 
balance  normal  to  the  oscillation  plane. 


* Resonance  occurs  from  equation  (72)  when 

* "nO'T) 


and 
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Loads  transmitted  to  the  model  are  deduced  from  strain-gage 
measurements  on  the  upper  and  lower  beams  shown  in  Figure  13, 
which  schematically  illustrates  the  forward  balance  arrangement. 
This  moment  beam  design  was  selected  to  provide  maximum  strain 
(i.e.,  gage  response)  with  minimum  deflection  so  that  the  model 
is  rigidly  forced  to  perform  the  motion  generated  by  the  shaft 
crank  and  yoke.  The  "degree  to  'which  this  is  achieved  is  related 
to  the  ratio  of  oscillation  frequency  to  the  mechanical  natural 
frequency  of  the  moment  beam  and  mass  located  forward  of  the 
beam.  Uote  will  be  made  of  this  point  when  balance  calibration 
and  data  reduction  are  discussed.  It  is  also  important  to  note 
that  the  moment  beam  is  located  between  the  pivot  and  the  model, 
and  therefore  loads  measured  by  the  gages  located  thereon  do  not 
contain  drive-train  or  pivot  frictional  components. 

The  instantaneous  angular  displacement  of  the  model  relative 
to  the  sting  is  obtained  from  strain  gages  mounted  on  the  center 
deflection  beam  mounted  between  the  model  and  the  sting,  as 
illustrated  in  Figure  13.  The  stiffness  of  this  beam  and  the 
inertia  of  the  model  and  pivoting  head  combine  to  generate  a 
natural  oscillation  frequency  of  the  balance  and  model.  This 
frequency  is  changed  by  the  addition  of  aerodynamic  stiffness 
during  wind-tunnel  flow  conditions,  but  typically  the  deflection 
beam  stiffness  is  greater  than  the  aerodynamic  stiffness,  as 
operation  at  as  high  a value  of  reduced  frequency  as  possible  is 
usually  desirable.  This  maximum  frequency  is  limited  by  the 
minimum  model  moment  of  inertia  ’which  can  be  practically  obtained 
and  by  the  maximum  frequency  allowable  by  both  the  data- 
acquisition  system  and  by  the  dynamic  characteristics  of  the 
balance  itself.  A foil  strain-gage  bridge  is  mounted  or.  the 
defelcticn  beam,  whereas  semiconductor  strain  gages  are  mounted 
on  the  moment  beams  to  provide  increased  sensitivity. 


In  order  to  obtain  satisfactory  performance  under  high- 
temperature  hypersonic  tunnel  flow  conditions,  the  balance  is 
provided  with  water  cooling.  However,  water  cooling  passages  are 
not  provided  within  the  strain-gage  beams  due  to  size  limitations, 
and  these  elements  can  be  subjected  to  some  temperature  draft 
during  testing.  The  bridges  are  temperature  compensated  and,  in 
addition,  since  oscillation  signals  are  the  primary  measurement 
objective,  a slow  bridge-zero  drift  due  to  uneven  thermal  loading 
can  be  eliminated  by  monitoring  only  the  a.c.  component  of  the 
bridge  outputs.  This  technique  does  not  account  for  bridge 
sensitivity  changes  due  to  uniform  temperature  changes;  however, 
account  can  be  taken  of  the  latter  by  calibration  and  by  appropriate 
design  tech.oiaue. 
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bridges  are  provided  on  the  beam  elements  to  develop  a normal- 
force  signal,  which  is  subsequently  combined  with  the  pitching 
moment,  analogous  to  a simple  moment- transfer  equation,  to  yield 
an  output  moment  signal  referenced  about  the  pivot  point.  This 
signal  combination  is  accomplished  on  the  balance  itself. 

Calibration  of  the  balance  and  associated  read-out  instru- 
mentation requires  consideration  of  both  static  and  dynamic 
effects.  Static  calibration  w as  accomplished  in  a conventional 
manner  by  dead-weight  loading  with  d.c.  bridge  excitation  voltages 
of  2,  *4,  5,  and  7 volts.  Actual  operation  of  the  balance 
was  performed  with  a bridge  voltage  of  aoproximately  li  volts. 


It  is  important  to  note  that  the  pitching-moment  beam 
deflections  under  actual  oscillating  conditions  represent  dynamic 
deflections  and,  as  such,  require  smaller  loads  for  equal  strains 
than  under  equivalent  static  conditions.  That  is,  the  instan- 
taneous beam  deflection  is  a function  of  applied  load  and  frequency, 
ana  in  principle  a dynamic  calibration  at  several  oscillation 
frequencies  is  required.  In  practice,  one  normally  assumes  that 
the  mechanical  natural  frequency  of  the  moment  beam  and  model 
inertia  combination  is  very  much  higher  than  the  oscillation 
frequency  during  testing,  and  hence  the  dynamic  deflections  can 
be  expected  to  be  approximately  equal  to  static  deflections  under 
equal  loading. 

An  indication  of  the  relative  dynamic  amplification  has  been 
obtained  by  earlier  tests  performed  on  the  balance  for  NASA.  A 
sketch  of  the  test  setup  for  that  study  is  shown  in  Figure  1^4. 
Sinusoidal  loads  were  generated  by  an  electro-mechanical  shaker 
and  applied  to  the  balance  through  a calibration  fixture.  The 
fixture  served  as  an  inertia  load  and  as  a means  for  load  application. 
A load  amplitude  cf  two  poinds  was  maintained  while  the  frequency 
was  varied  from  1 to  150  hertz.  The  balance  strain-gage 
bridge  outputs  were  recorded  and  analyzed  to  provide  the  frequency 
response  characteristics  of  the  balance-fixture  system.  The 
output  from  each  bridge  was  normalized  to  the  one  hertz 
reference  value. 

Figure  15  shows  a representative  frequency  response  curve 
for  the  pitching-moment  bridge  output  as  a function  of  load 
frequency.  A resonant  frequency  of  about  ^5  hertz  can  be  observed 
with  the  dynamic  response  ampli fication  being  about  10  at  this 
frequency.  .Mote  that  at  a frequency  of  aoout  10  hertz,  the  out- 
put is  approximately  10  percent  higher  than  the  static  value. 
Therefore,  during  actual  testing  at  this  frequency,  one  would 
reduce  the  output  measurement  about  10  percent  in  order  to  obtain 
the  actual  load  magnitude  from  static  calibration  information. 

This  result  is,  of  course,  dependent  on  the  fixture  inertia 
characteristics  being  representative  of  the  actual  model 
characteristics. 
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FIG.  14  SCHEMATIC  DIAGRAM  OF  BALANCE  FREQUENCY  RESPONSE  TEST  FIXTURE 
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Tiie  basic  data- reduction  analysis  for  forced-osci llation 
systems  has  been  described  and  is  not  repeated  here.  In  principle, 
if  measurements  are  made  of  the  applied  moment  amplitude, 
the  angular  displacement  amplitude ,«<  0 * the  phase  angle, 
between  the  moment  and  displacement  signals,  and  the  frequency  of 
oscillation,  cu,  then  the  total  system  damping,  C,  can  be  computed 
from  the  relation  (equation  (89)) 


c_  MoSinV 


u>«L 


(1?3) 


V.’hen  <j  is  approximately  equal  to  the  undamped  natural  frequency  of 
the  system,  the  phase  difference  becomes  approximately  90  degrees. 
At  this  frequency,  the  applied  moment  necessary  to  sustain  the 
motion  is  near  a minimum  value.  The  required  moment  is 
actually  a minimum  at  the  damped  natural  frequency,  which  is  slightly 
less  than  the  undamped  natural  frequency;  however,  the  difference 
between  the  two  is  negligible  for  lightly  damped  systems.  Wind- 
off  tare  measurements  were  made  under  near-vacuum  conditions 
prior  to  each  run  to  assess  the  structural  damping  influence. 

These  measurements  were  obtained  at  the  wind-off  resonant  fre- 
quency, which  was  typically  only  one  or  two  hertz  less  than  the 
wind-on  value,  due  to  the  large  deflee  i spring  stiffness 
relative  to  the  aerodynamic  stiffness. 


THE  WIND  TUNNEL  TESTS 


Tests  were  conducted  in  the  wind  tunnel  facilities  of 
Surface  Weapons  Cencer. 


Naval 


The  Hypervelocity  Research  Tunnel  ( 8 A ) operates  (Mach  18)  at 
temperatures  up  to  2400  Kelvin  in  the  supply  pressure  range  from 
200  to  700  atmospheres.  The  tunnel  utilizes  high-pressure  stored 
nitrogen  as  the  ’working  gas.  The  gas  is  heated  by  a forced- 
convection,  electrical  resistance-type  graphite  heater.  Under 
maximum  operating  conditions,  this  heater  delivers  power  to  the  gas 
at  a rate  of  1600  kilowatts. 


The  facility  provides  a maximum  Reynolds  number  of 
2 x 10s  per  meter,  which  represents  altitude  simulation  at 
36,000  meters.  The  average  steady-state  operating  time  at 
any  supply  condition  is  approximately  three  minutes. 

The  facility  is  supplied  with  a test  model  support 
mechanism  that  provides  pitch  and  roll  motions  and  readouts. 
Also,  instrumentation  and  data  processing  equipment  is 
available  for  force,  pressure,  temperature,  and  dynamic 
stability  testing. 
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Nozzle  flow  studies;  namely,  Pitot  pressure  surveys,  have 
been  made  at  supply  pressures  of  ^00  and  600  atmospheres. 

Results  of  these  studies  indicate  that  the  uniform  flow  test 
core  is  13  to  20  centimeters  in  diameter.  The  centerline 
Pitot  variation  over  a 71-centimeter  axial  length  is  approximately 
3 percent.  This  corresponds  to  a Mach  number  gradient  dM/dX  of 
about  0.5  x 10" 3 /centimeter. 

For  further  information  regarding  the  Mach  18  Hypervelocity 
Research  Tunnel,  see  MOLTR  7*  1—158 , Ref.  10. 10 

The  model  was  injected  into  the  Mach  18  flow  field  at  sting 
angle-of-attack  of  0°,  5°  and  10°,  disturbed  to  an  initial  angle-of- 
attack  and  permitted  to  freely  oscillate.  The  subsequent  induced 
strain  gage  signals  were  stored  on  magnetic  tape  and  converted  to 
digital  values. 

The  least  squares-di fferential  correction  technique  w as  then 
used  to  extract  the  stability  parameters  from  the  oscillatory 
motion. 


The  fitting  model  used  was 

oi-aKe^cosCtot**)*  d 


(12") 


The  stability  coefficients  determined  were  Cma,  the  damping 
moment  coefficient,  the  restoring  moment  coefficient,  and 
Cpj*  the  slope  of  the  restoring  moment  coefficient  with  angle-of- 


afctacx. 


[Tie  damping  moment  coefficient  was  determined  as: 


r _ £yi' 


The  Ditching  moment  was  determined  uS  • 


(125) 


C = 1±L_ 

* <*&«s 


The  slope  of  the  pitching  moment  coefficient  was 
determined  as: 


0 26) 


(127) 


5 Cornett,  R.  H.  Keel,  A.  G.,  "NOL  Mach  18  Hypervelocity 
Research  Tunnel, M Naval  Ordnance  Laboratory  Technical  P»epont 
7^-158,  August  197^ 
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where  the  assumption  has  been  made  that  £ and  £0  are  small  compared 
to  */TTl  » see  equation  (60). 

Result.-  from  a pitch  damping  investigation  conducted  in  Supersonic 
Tunnel  So.  ? and  the  Hypersonic  Tunnel  (8)  of  the  Naval 

weapons  Center,  White  Oak  Laboratory  are  also  presented 
:vi  :e  a Mach  number,  Reynolds  number  variation.  The 
rv/R:.  = .9",  10°  cone  was  run  at  a Mach  number  of  5 at  a free 
stream  Reynolds  number  of  9 x IQ6  per  foot.  For  a discussion  of 
the  operating  conditions  of  Supersonic  Tunnel  No.  2,  see 
Reference  11. 11  Tests  were  also  carried  out  in  the  Hypersonic 
Tunnel  (r)  at  a variety  of  stagnation  pressure  conditions.  For 
a discussion  of  the  operating  conditions  of  the  Hypersonic  Tunnel 
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A Mach  18  free  oscillation  data  point  was  also  obtained  in 
the  Hyper velocity  Research  Tunnel  for  a seven  percent  blunt 
a on  f i srura  t i on . 


WIND-TUNNEL  RES U LTS 

The  primary  objective  of  the  test  program  was  to  determine 
the  aerodynamic  stability  coefficients  for  a ten  degree  cone  as 
a function  of  bluntness  and  angle-of-attack  under  the  imposed 
testing  conditions,  and  to  examine  the  repeatability  of  these 
coefficients  with  angle-of-attack.  The  results  thus  obtained  are 
compared  with  static  test  results  as  well  as  results  from  other 
facilities  for  the  purpose  of  evaluating  the  data. 
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free  oscillation  system  possesses  sufficient  nensiti 
de  dynamic  stability  measurements  at  the  Hy serve loci 
Tunnel’s  operating  conditions.  The  reduceable  damp 
lens  thus  obtained  had  an  average  mean  amplitude  of 
1 cn  of  approximately  two  degrees  about  a trim  angle- 
Dti.ng  angle-of-attack  of  0°,  5°»  and  10°  were  used 
ut  the  testing  program.  The  flexure  pivot  with  the 
1 stiffness  was  selected  for  a majority  of  wind  tunn 
order  to  achieve  a larger  amplitude  of  oscillation. 
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Aero  and  Hydroballistics  Research  Facilities,"  Naval 


d.oance  Laboratory  Report  126 -'J , July  1967 


"^Geineder,  F. , Schlesinger,  M.  I.,  Baum,  G. , Cornett,  R., 
"The  U.  .G.  Naval  Ordnance  Laboratory  Hypersonic  Tunnel," 
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torsional  stiffness  provided  much  better  Derformance  in  that  very 
small  trim  angles  developed.  It  is  concluded  that  preliminary 
estimates  of  pitching  moment  should  be  made  prior  to  flexure 
selection  and  that  these  values  should  then  be  used  to  determine  the 
trim  angle  and  therefore  the  flexure  which  is  most  compatible 
with  the  test  requirements. 

No  frequency  effect  was  found  to  be  present.  The  system  was 
found  to  produce  highly  repeatable  data  yielding  a maximum 
deviation  of  3 percent  for  pitching  moment  coefficient  and 
8 percent  for  damping  moment  coefficient. 

Data  Reduction 

The  techniques  previously  described  were  used  to  fit  the 
aercballistic  theory  to  the  1-D  oscillation  obtained  from  the 
wind  tunnel  tests. 

The  free-oscillation  data  was  digitized  and  fit  each  record 
consisting  of  approximately  6000  data  points.  The  average  percent 
error  of  fit  of  the  theory  to  the  data  indicated  that  the  damped 
sinusoidal  equation.  Equation  (121),  represented  the  recorded 
wind  tunnel  motions  to  within  an  accuracy  of  2.5  percent.  A 
representative  plot  of  the  1-D  wind  tunnel  motion  is  presented 
in  Figure  16. 

Stability  Coefficients 

Representative  plots  of  pitching  moment  coefficient,  as 
determined  using  the  free  oscillation  technique,  versus 
angle-of-attack  are  given  in  Figures  17-19.  These  data  are  compared 
to  some  very  accurate  and  precise  force  and  moment  data13  taken 
for  the  same  model  in  the  Kyperveiocity  Research  Tunnel.  The 
results  of  the  force  and  moment  test  have  been  carefully  analyzed 
and  numerous  comparisons  have  been  made  with  other  experimental 
and  analytical  results.  Excellent  agreement  is  noted  between  the 
static  and  dynamic  results  over  the  bluntness  and  angle-of-attack 
range  considered.  Increase  in  angle-of-attack  yields  an  increase 
in  pitching  moment.  As  the  bluntness  is  increased  the  pitching 
moment  is  seen  to  increase  over  the  range  of  bluntness  considered. 


i "5 

"JKee 1,  A.,  "A  Note  on  the  Viscous  Effects  on  Center-of-Pressure 
Location  for  Sharp  and  31unted  Cones,”  Submitted  to  Journal  of 
Spacecraft  and  Rockets,  April  1975 
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FIG.  17  PITCHING  MOMENT  COEFFICIENT  VS  ANGLE-OFATTACK 
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FIG.  19  PITCHING  MOMENT  COEFFICIENT  VS  ANGLE  OF-ATTACK 
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In  order  to  further  evaluate  the  free-osei  nation  system 
comparison  of  Mach  18  free  oscillation  results  were  made  with 
free  oscillation  results  obtained  in  the  Aerospace  Research 
Laboratory  Mach  Facility. 1 “*  Results  typical  of  this  comparison 
are  given  in  Figures  20-25.  The  slope  of  the  pitching  moment 
curve  and  the  dynamic  stability  coefficient  Cm4! 

are  compared.  The  trends  exhibited  in  the  Mach  14  tests  are 
mirrored  by  the  Mach  13  results.  The  shifts  noted  between  the 
Mach  14  and  Mach  16  results  are  due  to  Mach  number  and 
small  center  of  gravity  location  effects.  The  effects  of  blunt- 
ness on  and  Cmi  are  presented  in  Figures  26  and  27, 

respectively.  The  center  of  gravity  location  referred  to  the 
virtual  nose  normalized  by  the  virtual  length  for  the  free 
oscillation  tests  had  a value  of  0.6078. 

FORGE  OSCILLATION  TEST  RESULTS 

Tests  were  carried  out  at  Mach  5 in  the  WOL  Supersonic  Tunnel 
(Tunnel  .*?)  and  Hypersonic  Tunnel  (Tunnel  "8)  at  a variety  of 
stagnation  pressure  conditions.  The  model  used  was  a 10-degree 
cone  having  a nose  bluntness  of  seven  percent.  The  analog  output 
signals  from  the  strain-gage  bridges  we re  digitally  sampled  and 
recorded  as  a function  of  time,  with  the  sampling  rate  being 
such  that  approximately  50  points  per  cycle  were  obtained.  On- 
line examination  of  the  pitching-moment  and  angular  displacement 
signals  displayed  on  a dual-trace  oscilloscope  provided  the 
necessary  feedback  information  for  confirming  resonant  conditions. 
Manual  motor  speed  control  proved  to  be  satisfactory  for  these 
tests.  The  phase  difference  between  signals  changed  dramatically 
as  motor  speed  increased  through  the  desired  resonant  condition, 
and  the  moment  signal  minimization  at  resonance  was  quite  definite. 
Post-test  examination  of  the  recorded  signals  indicated  that  a 
phase  difference  within  a few  degrees  of  the  90-degree  resonance 
value  was  normally  achieved.  Checks  were  also  made  to  insure 
that  phase  shifts  within  the  instrumentation  were  identical  for 
both  signals  of  interest. 

The  more  comprehensive  testing  was  accomplished  in  the  Super- 
sonic Tunnel,  with  data  being  obtained  at  pitch  angles  of  0,  fc,  and 
11  degrees;  yaw  angles  of  0 and  7 degrees;  and  at  total  pressures 
of  30,  30,  100,  150,  and  200  psia.  After  operation  of  the  balance 


“ ’V'alchner,  0.,  Sawyer,  F.  M. , "In-Plane  and  Out  of  Plane 
Stability  Derivatives  of  Slender  Cones  at  Mach  14,"  Aerospace 
Research  Laboratories  Report  73-0090,  July  1973 
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and  data-acquisition  system  was  confirmed  in  the  Supersonic  Tunnel, 
a brief  test  was  also  conducted  in  the  Hypersonic  Tunnel  at 
Mach  5 and  at  a total  pressure  of  370  psia  to  validate  the 
balance  and  support  system  performance.  Some  of  the  results  of 
these  tests  are  presented  in  Figures  28-31  and  compared  with 
available  data. 


Figure  28  shows  the  measured  damping  coefficients  as  a function 
of  Reynolds  number  and  compares  these  results  with  some  ballistic 
range  data  at  Mach  5-7  to  7.5  (Ref.  11)  and  with  some  free- 
oscillation  wind-tunnel  data  obtained  at  Mach  6. 15  There 
appears  to  be  little  Reynolds  number  influence  on  the  data 
presented.  The  lowest  Reynolds  number  at  which  the  present  data 
•were  obtained  is  representative  of  the  minimum  value  at  which  the 
damping  moment  could  be  accurately  measured  at  a bridge  excitation 
voltage  of  four  volts  and  with  the  present  instrumentation  setup. 
Balance  sensitivity  to  lower  loads  can  be  achieved  by  increasing  the 
bridge  excitation  (this  is  limited  by  the  requirement  for  minimizing 
Joule  heating  of  the  gages  and  by  voltage  limitation  of  the  gages 
themselves)  and  by  increasing  the  signal  gain  and  filtering 
characteristics  of  the  data-acquisition  and  recording  equipment. 

Figure  29  presents  a similar  comparison  of  damping  data  as  a 
function  of  Mach  number.  The  figure  includes  a Mach  18  data  point 
obtained  in  the  'a'OL  Hypervelocity  Research  Tunnel  by  the  small- 
amplitude  free-osclllation  technique  for  the  7 percent  blunt  case. 
Both  the  Mach  5 data  and  Mach  13  data  appear  to  be  in  reasonable 
agreement  with  available  data.  It  should  be  noted  that  several 
adjustments  could  be  made  to  the  data  to  account  for  center-of- 
osci llatior.  and  nose  bluntness  difference  but  these  corrections  are 
relatively  small  and  are  not  warranted  for  the  discussion  here. 


The  in-plane  dampign  data  are  presented  in  Figure  30  as  a 
function  of  mean  angle  of  attack  and  compared  with  some  Mach  2 
10-degree  cone  (R  /R  =0.017)  data.  The  Mach  2 data  were  obtained 
with  a single-degree-of-freedom,  free-oscillation  balance 
incorporating  a cross-flexure  pivot.16  The  present  Mach  5 data 
and  the  Mach  2 data  indicate  no  significant  angle-of-attack 
effects  for  angles  less  than  the  cone  half  angle.  An  out-of-plane 
damping  data  point  obtained  at  an  angle  of  attack  of  seven  degrees 
with  the  present  forced-oscillation  technique  (C.**  + 
id) 


;s  obtained)  also  agrees  in  a similar  manner  with  the  Mach  2 


--  °-»;1ane  damning  data  of  Reference  lb.  However,  the  out-of- 
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damping  at  an  angle  of  attack  near  the  cone  half  angle,  but 
rather  they  indicate  some  decrease  in  damping  for  this 
configuration  as  angle  of  attack  approaches  30  degrees  or  so. 


Finally,  a data  point  was  obtained  at  zero  mean  angle  of 
attack  in  the  ’.v'L  Hypersonic  Tunnel  for  the  7 percent  blunt  case. 

The  results  obtained  were  in  agreement  with  the  data  obtained  in 
the  Supersonic  Tunnel.  The  balance  support  system  in  the  Hypersonic 
Tunnel  dees  not  permit  out-of-plane  damping  data  to  be  obtained 
with  this  particular  balance;  however,  in-plane  damping  may  be 
obtained  at  angles  of  attack  up  to  30  degrees , a limit  imposed  by 
the  necessity  for  maintaining  an  adequate  oil  both  for  the  motor 
gearing  in  the  drive-unit  housing  of  the  balance,  discussed 
earlier.  The  allowable  load  limits  imposed  by  the  mechanical 
design  of  the  balance  were  not  reached  in  any  of  the  tests 
reported  here,  and  therefore  the  balance  possesses  the  capability 
for  use  at  higher  loads  (i.e.,  higher  total  pressures)  than 
indicated  by  these  evaluation  tests. 


RANGE  WIND-TUNNEL  PITCH  DAMPING  DISCREPANCIES 


Reference  1 reports  an  increase  in  dynamic  stability  as 
Mach  number  is  increased.  As  a result  of  this  range  conducted 
investigation,  it  was  concluded  that  the  measurements  made  suggest 
that  the  damping  derivatives  at  higher  Mach  numbers  may  have 
some  dependence  on  the  type  of  motion  pattern  experienced  by  the 
model.  As  previously  seen,  in  comparing  the  data  of  Reference  1 
with  comparable  wind  tunnel  data,  the  wind  tunnel  determined 
values  indicate  a decrease  in  dynamic  stability  with  Mach  number, 
which  is  not  in  agreement  with  the  results  of  the  range  study. 

No  Reynolds  number  variation  can  be  seen  to  contribute  as  was 
shown  in  Figure  23.  Reference  1 also  concludes  that  the  reported 
discrepancies  may  be  due  to  model  support  interference.  During 
the  present  investigation,  no  conclusive  evidence  of  this  type  of 
interference  could  be  determined.  The  high  repeatability  of  the 
data  and  the  excellent  agreement  noted  with  the  static  test 
results  and  results  from  other  facilities  4 end  to  de-emphasize 
this  conclusion.  It  is  felt,  however,  that  this  discrepancy  may 
be  explained  by  a difference  between  the  in  plane  and  out  of 
plane  dynamic  stability  coef fi cients . Such  results  were  lirst 
rooorted  for  angles-of-attack  greater  than  the  cone  half  angle 
(Ref.  16). 

Tabok  and  Scniff17  suggested  a theoretical  moment  formulation 
for  nonplanar  vehicle  motion.  This  formulation  was  nonlinear  ’with 
respect  to  angular  displacement  and  linear  with  angular  rates. 


1 ^Tohak , M.,  Schiff,  L. , "A  Non line 
Formulation  and  Its  Implication  for  Dyn 
NASA  Ames  Research  Censer,  AIAA  Paper  7 
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It  was  concluded  that  yaw  damping  at  angle-of-attack  must  be 
different  from  pitch  damping  if  the  pitch  moment  is  a nonlinear 
function  of  angle-of-attack.  It  has  been  shown  that  slightly 
blunted  cones  have  highly  nonlinear  aerodynamics  at  small 
angles-of-attack , which  are  only  a fraction  of  the  cone  half 
anrle  (References  14  and  IS)18  and  subsequent  analysis  followed.19 

Two- degree- of- freedom  tests  have  also  been  conducted20  to 
investigate  the  in-plane,  out-of-plane  defferenee.  It  was  found 
that  a transverse  product  of  inertia  has  a very  strong  coupling 
effect  in  the  case  of  nonplanar  motion.  In  the  case  of  planar 
motion  linear  aerodynamics  were  experienced.  In  the  case  of 
nonplanar  motion  an  aerodynamic  accleration  coupling  derivative 
was  extracted  which  had  the  same  effect  as  an  inertial  coupling 
by-product  of  inertia.21 


Such  damping  asymmetry  tends  to  invalidate  the  use  of  the 
tri-cyclic  theory  to  correctly  represent  the  free  flight  motion 
of  blunted  sphere  cones.  Out-of-plane  damping  is  then  present 
in  the  range  test  and  is  not  present  in  the  wind  tests. 
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If  the  range  data  were  fit  without  allowing  for  this  in-plane,  out- 
of-plane  difference  an  error  will  be  introduced  which  could  lead 
to  the  over  prediction  of  total  damping.  As  the  Mach  number  is 
increased  the  initial  angular  disturbance  of  the  range  model  would 
increase  and  thus  the  resulting  oscillations  and  the  over  prediction 
of  total  dynamic  stability  coefficient. 

The  reported  in-plane,  out-of-plane  damping  difference  should 
be  further  investigated.  A corrected  equation  of  motion  should 
be  determined  *nd  simulation  conducted  to  determine  the  overall 
effect  of  such  coefficient  asymmetry.  In  the  meantime,  care 
should  be  taken  in  comparing  results  determined  from  methods 
where  in-plane,  out-of-plane  damping  is  present  and  methods 
where  only  ir,-piane  damping  is  present  until  the  source  for  the 
disagreement  is  well  understood. 

SUMMARY 


The  damping  in  pitch  characteristics  of  a ten  degree  cone 
have  been  investigated.  Numerous  single-degree-of- freedom  wind 
tunnel  experiments  have  been  conducted  in  order  to  obtain  the 
variation  of  damping  coefficient  with  bluntness,  Mach  number, 
Reynolds  number  and  angle-of-attack . The  results  of  the  present 
investigation  have  been  compared  with  results  from  other 
facilities  and  have  been  found  to  be  in  good  agreement.  Data 
on  the  ten  degree  Supersonic  Tunnel  Association  cone  was  also 
obtained  at  Mach  18.  Damping  data  on  this  configuration  at 
Mach  13  was  not  previously  available.  A discrepancy  was  noted 
with  damping  data  variation  with  Mach  number  as  obtained  from  a 
ballistic  range  experiment.  An  explanation  of  this  discrepancy, 
based  on  in-plane,  out-of-plane  damping  difference,  has  been 
offered.  A detailed  explanation  of  the  theory,  data  reduction 
technique  and  experimental  methods  has  also  been  presented  in 
order  to  document  the  dynamic  testing  capabilities  of  the  Naval 
Surface  Weapons  Center,  White  Oak  Laboratory. 


? 


